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Within the impulse approximation the response of a many body system at large momentum tranfer 
can be written in a simple and transparent form, allowing to directly relate the inclusive scatter- 
ing cross section to the properties of the target ground state. Although the physics assumptions 
underlying impulse approximation are well defined, their implementation involves ambiguities that 
may cause significant differences in the calculated responses. We show that, while minimal use of 
the impulse approximation assumptions naturally leads to write the response in terms of the target 
spectral function, the widely used alternative definition in terms of the momentum distribution in- 
volves a more extended use of the same assumtpions. The difference between the responses resulting 
from the two procedures is illustrated by two examples. 
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In this letter we show that a truly unambiguous form of 
the response of a strongly interacting many-body system 
to an external probe, within the impulse approximation, 
is based on the use of the spectral function, rather than 
the momentum distribution. 

The main assumption underlying impulse approxima- 
tion (lA) is that, as the spacial resolution of a probe de- 
livering momentum q to a many body system is l/|q|, 
at large enough |q| the target is seen by the probe as a 
collection of individual constituents. Within this picture, 
the response measures the probability that, after giving 
one of the constituents a momentum q at time t = 0, 
the system be reverted to the ground state after time t 
giving the same constituent a momentum — q. 

The second assumption involved in lA is that final 
state interactions (FSI) between the hit constituent and 
the (N-l)-particle spectator system be negligible. The 
most popular argument supporting this assumption is 
based on the observation that, compared to the ampli- 
tude in absence of FSI, the amplitude of the process in- 
cluding a rescattering in the final state involves an extra 
propagator, describing the motion of the struck particle 
carrying a momentum ~ q. As a consequence, this am- 
plitude is expected to be suppressed when |q| is large. 

In spite of the fact that the two basic assumptions un- 
derlying lA can be unambiguously stated, in the litera- 
ture one finds two different definitions of the I A response, 
involving either the target spectral function or its mo- 
mentum distribution ||^. 

The two different definitions arise from different im- 
plementations of the I A assumptions, and may lead to 
significantly different numerical results. In addition, as 
lA can be seen as the zero-th order of a systematic ap- 
proximation scheme, to be improved upon including FSI 
effects, the ambiguity in the lA response poses a serious 



problem, making it difficult to identify genuine FSI ef- 
fects. This feature is particularly critical in the analysis 
of the electromagnetic response of nuclear systems, where 
FSI are believed to play a relevant role even at large |q| 

i- 

In this short note we show that the definition of the re- 
sponse in terms of the spectral function follows from min- 
imal use of the assumptions involved in the lA scheme, 
and correctly takes into account the correlation between 
momentum and removal energy of the hit constituent. 
On the other hand, a more extended use of the same 
assumptions leads to the definition in terms of the mo- 
mentum distribution, that totally disregards the effect of 
the removal energy distribution. 

The response of a N-particle system to a scalar probe 
is defined as 

5(q,^) = i J ^ e-*(0|pt(t)p^(0)|0) 

= k I ^ e^^*(0|e«V^e-«Vq|0) , (1) 

where q and lo denote the momentum and energy trans- 
fer, respectively, H and |0) are the target hamiltonian 
and ground state, satisfying the Schrodinger equation 
HIO) = Eo|0), and pq = J2k"'k+q^^^^ «k+q and aw be- 
ing the usual creation and annihilation operators. Note 
that the above definition can be readily generalized to 
describe the electromagnetic response replacing pq with 
the appropriate electromagnetic current operator. 

Using the Schrodinger equation to get rid of one of the 
propagators appearing in eq.(|^) we obtain 

5(q,c.) = i J |- e^(-+^°)M0K e-HVq|0) . (2) 

The above definition can be simplified introducing the 
first assumption involved in lA, i.e. that the process in- 
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volves only one constituent, while the remaining (N-1) 
particles act as spectators. As a result, the ground state 
expectation value appearing in eq.(^ can be rewritten 
in configuration space as (i? = (ri, . . . , tn) specifies the 
positions of the N target constituents) 

{0\pI e-*"Vq|0) =N J rfi?di?'*5(i?)e-*i"-i 

x(i?|e~*"*|i?>'i"-'i^'o(i?/) , (3) 

where ^o(-R) — {R\0) is the ground state wave function. 
The N-particle hamiltonian H can be split according 

to 



H = Ho + Ti + H 



FSI 



(4) 



where Hq denotes the hamiltonian of the spectator sys- 
tem 
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(5) 
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Vij and m being the potential describing the interactions 
between target constituents and the constituent mass, 
respectively. The remaining two terms in eq.(^ are the 
kinetic energy of the struck particle. 



V? 

2m ' 



and 



N 



Hfsi = E • 

J=2 



(6) 



(7) 



The second assumtpion involved in lA amounts to dis- 
regard the contribution of Hpsi, describing the FSI be- 
tween the hit constituent and the spectators. As Hg and 
Ti obviously commute, this allows to rewrite the config- 
uration space N-body propagator appearing in eq.(||) in 
the simple factorized form (i? = (r2 . . . tn)) 



(i?|e-'"*|i?') = (i?|e-*"°*|i?')(i-i|e""'*|ri) 



(8) 



The two propagators in the rhs of the above equation 
can be written in spectral representation as 



{R\e-'^"'\R') = ^e-*^"*$„(i?)$;(i?') 



and 



-iTit|„/ 



(2^ 



-iBpt ip-(ri-r'i) 



(9) 



(10) 



where $„(-R) = {R\n), En and |n) satisfy the (N-1) 
particle Schrodinger equation Ho|n) = i?„|n), and Ep - 
/2m. 



Using eqs.(^) and (^0|) and substituting eq.(^) into 
eq.(||) we get 



-i{Ep+E„)t 



dR e*(P-i) ''i^'*(i?)$„(^) 



(11) 



Finally, substitution of the above result into eq.(y) 
leads to 

5(q,c.)=. J ^dE P{p-ci,E) 

xS{lj ~ Ep~ E) , (12) 
where the spectral function, defined as 

P(k,E) e* '-i«'S(i?)<i>„(i?) 

xS{E + Eo~En), (13) 

measures the probability of removing a constituent of 
momentum k from the target ground state leaving the 
residual system with excitation energy E. 

Let us now consider a different way of implementing 
the physical assumptions underlying lA in the calcula- 
tion of S{q,uj). In going from eq.(|^) to eq.(§) we have 
exploited Schrodinger equation to get rid of one of the two 
N-body propagators. We have then used the assumption 
Hp SI = to rewrite the remainig propagator in the fac- 
torized form that led to the emergence of the spectral 
function in the formalism. In principle, since lA provides 
a prescription to rewrite the N-particle propagator in a 
simpler form, one may just as well use this prescription 
and rewrite both propagators appearing in eq. (H) , rather 
than use Schrodinger equation. However, this procedure 
results in a definition of >S'(q, lu) in which the information 
on the target removal energy distribution is totally lost. 

The ground state expectation value relevant in this 
case. 



Pq|0) J dRdR'dR"%{R) 
(i?|e™*|i?"}e~*'i-'""(i?"|e"™*|i?')e"'"''^*o(i?') , (14) 



can be rewritten using again factorization and the spec- 
tral representation. In addition, the dependence upon 
the state of the spectator system can be removed apply- 
ing the orthonormality relations 



and 



di?$;(i?)$„(i?) = s„ 



J2Kimn{R') = SiR-R') 



(15) 



(16) 
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As a result, the rhs of eq.(|lj) becomes 



n j dR d^r[ d^r'l j ■ 



e'(-^'=-^'')***(ri,i?) 



(27r)3 (27r)3 
and integration over r" and p yields 

X / dVidV'ie^ t'-i-'-'i) / **(ri,^)*o(r'i,^) • (18) 



Finally, substitution of the above equation into eq. (|^) 
leads to 



^(q,c.) 



(2^ 



7i{\^) 5{lo + Ek ~ E\^+^{) , (19) 



where the momentum distribution n(k), yielding the 
probability to find a constituent carrying momentum k 
in the target ground state, is given by 



i(k) dVidVie'*'('"i-'"'i) 

X [ dR *5(ri,i?)*o(r'i,i?) 



(20) 



Comparison between the above equation and eq-dl^) 
shows that the momentum distribution is simply related 
to the spectral function through 



n(k) = J dE P{k, E) 



(21) 



As a first example, illustrative of the differences be- 
tween ^(q, w) evaluated using eq.(^2|) and that resulting 
from eq. ([l9|) , we will discuss the response of infinite nu- 
clear matter at equilibrium density p — 0.16 fm"^. 

An ah initio microscopic calculation of the nuclear mat- 
ter spectral function, carried out within the framework 
of Correlated Basis Function (CBF) perturbation theory 
using a realistic nuclear hamiltonian, is described in ref. 
10 . The main feature of the spectral function of ref. ||l| is 
the presence of a substantial amount of strength at large 
E, leading to an average removal energy e = {E) = 61.9 
MeV, much larger than the Fermi energy = 16 MeV. 
In addition, the calculated P(k, E) exhibits a strong cor- 
relation between momentum and removal energy, imply- 
ing that large momentum (|k| >> kp, kp — 1.33 fm^^ 
being the Fermi momentum) always corresponds to large 
removal energy [E » ep). For example, 50% of the 
strength at |k| = 3 fm^^ resides at > 200 MeV 0. 

The solid and dashed lines of fig.l show the uj depen- 
dence of S'(q, w) evaluated from eqs.(|l2|) and (|l9|), re- 
spectively, at |q|= 5 fm^^. At this momentum transfer, 
the nuclear response exhibits scaling in the variable y , 
reflecting the onset of the lA regime |^ . 



The solid line of figJi| has been obtained using the spec- 
tral function of ref. ||ll]7 whereas the momentum distribu- 
tion employed to obtain the dashed line has been consis- 
tently calculated by E integration of the same P(k, i?), 
according to eq.([2l|). 

While the two curves have similar shape, their width 
being dictated by the momentum distribution, they ap- 
pear to be shifted with respect to one another. The peak 
of the dashed curve is located at energy w ~ |qp/2mAr, 
corresponding to elastic scattering off a free stationary 
nucleon, whereas the solid line, due to the removal energy 
distribution described by the spectral function, peaks 
at significanlty larger energy. To illustrate this feature 
we show by diamonds the results obtained shifting the 
dashed curve by e = 61.9MeT^. 

In addition to the shift in the position of the peak, 
the dashed and solid curves sizably differ at low energy 
tranfer, where the response obtained using the momen- 
tum distribution is much larger than that obtained from 
eq.(|lj). The difference between the two curves in the 
low to region points to the fact that corrections to the re- 
sponse obtained from eq.(|9|) cannot be unambiguously 
identified as arising from mechanisms not included in the 
definition of lA. For example, a quantitative study of FSI 
effects, that are known to dominate the nuclear response 
at low w, can only be carried out starting from S{q,uj) 
defined as in eq.(p^. 
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FIG. 1. Infinite nuclear matter 5'(|q|,aj) at equilibrium 
density and [q| = 5 fm~^. The solid and dashed lines have 
been obtained from eqs. (|l|) and (^, respectively. The di- 
amonds represent the results obtained shifting the dashed 
curve by e=61.9 MeV. 

The results of fig.0 clearly show that the nuclear re- 
sponses extracted from electron-nucleus scattering data 
at momentum transfers in the few GeV/c range (1 GeV/c 
^ 5 fm~^) must be analyzed using spectral functions, as 
in !§. 

On the other hand, the definition of ^(q, oj) in terms 
of the momentum distribution has been successfully used 
to describe the response of liquid helium, measured by 
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inclusive scattering of thermal neutrons y|7|]. The ex- 
cellent agreement between the response calculated from 
eq.(|l^) and the experimental one can be explained not- 
ing that i) the region of momentum transfer covered by 
neutron scattering data extends to extremely high |q|, 
typical values being larger than 10 A~^, and ii) the anal- 
ysis has been focused on the region of the peak. 

In Hquid '^He at equihbrium density (p = 0.01635 A^"^) 
the half-width of the peak of the response at |q| = 10 A^^ 
is roughly given by (M denotes the mass of the helium 
atom) Iql/cp/M ~ 200 °K, to be compared to a Fermi 
energy ep = 2.47 °K, and the shift in w of - 10 °K 
produced by the removal energy of the struck particle 
reduces to a very small effect. 

The nuclear matter response of figj^, on the other 
hand, has a half-width of ~ 250 MeV, to be compared to 
a Fermi energy ep = MeV and an average removal en- 
ergy e — 61.9 MeV. As a consequence, the shift between 
the solid and dashed lines is clearly visible j|] . To observe 
a comparable effect in liquid ^He, one should consider the 
response at |q| '^3 A^^, where the half- width of the peak 
shrinks to ~ 50 °K. 
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FIG. 2. 5'(|q|,a;) in liquid ^ He at |q| = 10 A"^ and equilib- 
rium density. The solid and dashed lines have been obtained 
using eqs.(|l2[) and (|l^), respectively. 

The small effect of the removal energy on the position 
of the peak of the response of liquid ^He at p = 0.01635 
A~'^ and |q| = 10 A~^ is illustrated in fig.^, where the 
solid and dashed lines correspond to S'(q, us) evaluated 
from eqs.(^2|) and (p^), respectively. The momentum 
distribution and spectral function employed in the calcu- 
lations have been consistently obtained within the Fermi 
Hypernetted Chain formalism and the Diffusion Monte 
Carlo method |^ . 

In conclusion, we have shown that the two different 
prescriptions used in the literature to evaluate the re- 
sponse of a strongly interacting many body system cor- 
respond to different implementations of the assumptions 
underlying lA. 

While the definition in terms of spectral function re- 



quires minimal use of the approximations and correctly 
takes into account the removal energy distribution of the 
struck particle, the response obtained from the momen- 
tum distribution does not include all interaction effects. 

As shown by the excellent agreement between theory 
and the experimental data for liquid '^He at |q| > 10 
A~^ Q, this feature does not appear to be critical to the 
analysis of the response at very large momentum transfer, 
corresponding to (|q|/fcF) > 10, in the region of the peak. 

On the other hand, away from the peak large discrep- 
ancies between the S'(q, w) obtained from eqs.(|l2|) and 
( p9| ) persist even at very large |q|. Hence, a quantitative 
study of FSI effects, that are known to be important in 
the low energy region uj << |q|^/2m, requires as start- 
ing point the lA response calculated using the spectral 
function. 
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